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ABSTRACT 

The accretion disk in AGN is expected to produce strong outflows, in particular 
a UV-line driven wind. Several observed spectral features, including the soft X-ray 
excess, have been associated with the accretion disk wind. However, current spectral 
models of the X-ray spectrum of AGN observed through an accretion disk wind, known 
to provide a good fit to the observed X-ray data, are ad-hoc in their treatment of the 
outflow velocity and density of the wind material. In order to address these limitations 
we adopt a numerical computational method that links a series of radiative transfer 
calculations, incorporating the effect of a global velocity field in a self-consistent man- 
ner (XSCORT). We present a series of example spectra from the XSCORT code that 
allow us to examine the shape of AGN X-ray spectra seen through a smooth wind 
with terminal velocity of 0.3c, as appropriate for a UV-line driven wind. We calculate 
spectra for a a range of different acceleration laws, density distributions, total column 
densities and ionization parameters, but all these have sharp features that contrast 
strongly with both the previous 'smeared absorption' models, and with the observed 
smoothness of the soft X-ray excess. This rules out absorption in a radiatively driven 
accretion disk wind as the origin of the soft X-ray excess, though a larger terminal 
velocity, possibly associated with material in a magnetically driven outflow /jet, may 
allow outflow models to recover a smooth excess. 

Key words: quasars: general radiative transfer accretion, accretion discs X-rays: 
galaxies galaxies: active galaxies: Seyfert 



1 WHY ARE WINDS IMPORTANT? 

Mass outflows from accretion flows are observed to be ubiq- 
uitous features in accretion-powered astrophysical sources, 
with winds (and jets) being seen in systems as diverse as 
AGN and young stellar objects (e.g. the review by Livio 
2003). These outflows are expected to carry both energy and 
angular momentum away from the associated accretion disk 
material and, as such, are intimately linked to the accretion 
flow (e.g. Blandford & Payne 1982, Ampleford et al. 2007, 
Namouni 2007), however, the driving mechanism, particu- 
larly in the case of winds, may vary considerably between 
different sources. Different driving mechanisms give different 
physical properties (mass loss rate, outflow velocity, angu- 
lar momentum) for the wind material, and thus the impact 
of the wind on the underlying accretion flow, as well as on 
the immediate environment, will be considerably different 
for different sources. 

The most plausible physical mechanisms for launching 
a wind from an accretion disk around a black hole are ther- 
mal driving, radiation driving (both line & continuum) and 



magnetic driving. Thermal driving results from irradiation 
of the outer disc by emission from the bright inner accre- 
tions flow and the coronal X-ray source. This creates a pres- 
sure difference between the irradiated surface of the disk 
and the vacuum above, which drives material up to larger 
scale heights. The Compton temperature is typically ~10 7-8 
K and is constant across the disk, however the escape ve- 
locity decreases with radius. Beyond a certain radius, then, 
the material driven out of the disk will have sufficient ther- 
mal energy to exceed the escape velocity of the system at 
that point. This launch radius, typically ~10 4 Schwarzschild 
radii, is given by Riaunch>W 10 (M BH /Tc.s) cm, where M B h 
is the mass of the black hole in solar masses and Tc,s is 
the Compton temperature in units of 10 8 K (Begelman et 
al. 1983). 

Thermal winds should occur in all discs that have suf- 
ficient mass beyond the launch radius and so represent the 
minimum wind that is expected to be present in accreting 
sources. The mass loss rate of these winds can be dramati- 
cally enhanced if this material above the disc absorbs a sub- 
stantial amount of radiation from the central source, with 
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the imparted outward momentum giving rise to rapid radial 
acceleration of the wind material. Accounting for the ra- 
diation pressure, the effective gravity binding the material 
to the black hole is then reduced by a factor (1-f L/L Bdd ). 
Here, L is the source luminosity, L_gdd is the Eddington lu- 
minosity and f—(a abs +aT)/u T where a abs and <jt are the 
true absorption and electron scattering cross-sections, re- 
spectively. For highly ionized material, f is dominated by 
the electron scattering cross-section which thus represents 
the minimum radiation pressure we expect to be exerted 
on any wind material, and gives rise to continuum radia- 
tively driven winds around L/Lsdd~l. For less highly ion- 
ized material, f becomes dominated by the true absorption 
cross-section resulting in absorption line driving which can 
accelerate a wind at considerably lower L/Lem- 

The potential of magnetic fields to drive a wind is dif- 
ficult to quantify because the properties of such a wind 
will depend strongly on the field configuration (Blandford 
& Payne 1982; Proga 2000; 2003). Some progress is being 
made towards understanding this {e.g. Hawley & Krolik 
2006, Kato 2007), but as yet there is no clear consensus, 
beyond the general agreement that winds will be substan- 
tially enhanced by magnetic fields, and that such a wind 
could be launched from any radius in the disk. This lack of 
diagnostic power means that magnetic driving can only be 
implied as the dominant wind driving mechanism for sources 
in which all the other mechanisms have been ruled out. 

Moderate Eddington fraction [i.e. ~0. 1-0.5) stellar 
mass black holes and neutron stars, seen at high inclinations, 
often show high ionization blue-shifted absorption lines, typ- 
ically associated with H-like and He-like iron. These lines 
imply the presence of columns of ifj 23 ~ 24 cm~ 2 , outflowing 
at ~200-500 km s _1 . Continuum radiation is ruled out as the 
dominant driving mechanism behind this material because 
the luminosities are considerably less than Eddington, while 
the high ionization state of the wind material rules out line 
driving. Instead the observed properties are generally well 
matched to thermal wind models {e.g. Woods et al. 1996), 
however there remains considerable debate in the commu- 
nity regarding the role of additional driving mechanisms, 
particularly concerning the role of magnetic fields (Miller et 
al. 2006, Netzer 2006). 

The situation is slightly different in AGN due to the 
lower characteristic temperature associated with an accre- 
tion disk around a larger mass black hole. Coronal X-rays 
again push the material in the disk towards ~10 7 K, implying 
that a thermal wind could be launched from a similar grav- 
itational radii in the disk as the thermal winds observed in 
X-ray binaries. However the lower disk temperature means 
that material lifted from the disk to larger scale heights, will 
not necessarily be as highly ionized as in the case of stellar 
mass black-holes, and that the accretion disk will produce 
a copious flux of UV photons. Such partially ionized ma- 
terial has huge opacity to UV and soft X-ray photons and 
resonance line absorption in such material can thus power a 
stronger wind, from smaller radii, than thermal driving. As 
such, UV line driving is expected to be the foremost physical 
mechanism for AGN accretion disk winds (Murray & Chiang 
1995, Proga et al. 2004). The predicted geometry for a UV- 
line driven wind has three distinct regions; a low density, 
completely ionized, funnel perpendicular to the disc, a slow, 
mostly neutral, outflow concentrated towards the plane of 



the disk, and a fast (~0.1-0.2c), partially ionized, outflow 
in-between, forming a bipolar funnel {e.g. Proga 2004). 

Observationally, there is clear evidence for winds on 
large scales in AGN. In the X-ray band, a wide range of 
atomic transitions from partially ionized material (the warm 
absorber systems) are observed in many type-1 AGN {e.g. 
Kaastra et al. 2002, Blustin et al. 2005), however the inferred 
radii of this material places it in the region of the putative 
molecular torus, considerably further out than winds origi- 
nating from an accretion disk, suggesting an origin for this 
material in a wind from an irradiated torus rather than from 
the accretion disc itself (Krolik & Kriss 2001). Similarly, at 
higher inclination angles, intermediate type AGN show clear 
evidence for large column densities of moderately ionized ab- 
sorbing material at large radii (most notably in NGC 4151 
and Mrk 6 - Schurch & Warwick 2002; Schurch, Griffiths 
& Warwick 2006), however it is, again, unclear whether this 
material can be linked with winds from an accretion disk. In 
type-2 AGN, signatures from the accretion disc, and features 
from any associated wind, are overwhelmed by absorption 
in the neutral material of the obscuring torus. 

There is also evidence for outflows with properties sim- 
ilar to those we expect from the putative UV-line driven 
accretion disk wind. BAL QSO's, for example, show highly 
ionized absorption lines in their optical and X-ray spectra, 
with relativistic outflow velocities of up to ~0.2c (Weymann 
et al. 1991, Pounds & Page 2006), and it has been suggested 
that the putative molecular torus may, in fact, be either 
the slowly outflowing, high column density, neutral material 
near the plane of the accretion disk, or an accretion disk 
wind that has stalled at large radii (Elvis 2000). 

More indirect evidence for an outflow from the accre- 
tion disk comes from observations of the continuum spectral 
shape in high L/L B(id AGN, which typically shows a soft X- 
ray excess at energies below ~1 keV {e.g. Porquet et al. 2004, 
Crummy et al. 2006). This smooth component shows no ob- 
vious atomic features, and can be modelled with a thermal 
continuum. Fitting a thermal continuum to a range of AGN 
X-ray spectra reveals that the soft excess has a characteris- 
tic 'temperature' that is remarkably constant across a large 
number of sources, despite the large range of intrinsic source 
properties, making it unlikely that the soft X-ray excess is 
a true continuum component. The simplest mechanism for 
producing such a fixed 'temperature' component is atomic 
transitions in partially ionized material, in particular those 
associate with the strong jump in opacity between ~0.7-3 
keV caused by O VII/ VIII and absorption in partially ion- 
ized iron. Fabian et al. (2002) suggested that the soft excess 
could be linked to partially ionized material in the optically 
thick accretion disc seen in reflection, with the relativistic 
motion of the disc smearing the characteristic sharp atomic 
features into a smooth, continuum-like, component. Alterna- 
tively, Gierlinski & Done (2004 - hereafter GD04) suggested 
that the soft excess arises from absorption (and associated 
emission: Schurch & Done 2006 - hereafter SD06) in partially 
ionized, optically thin material. In this case, the large veloc- 
ity shear required to smooth the otherwise sharp atomic fea- 
tures implies that this material forms a wind and/or highly 
turbulent structure above the disc. 

Remarkably, this picture suggests that the accretion 
disk wind may be the dominant force shaping the X-ray 
spectrum we observe, and in some cases (notably 1H 0707- 
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495; Sobolcwska & Done 2007), the impact of the wind may 
be so severe that it completely masks the underlying con- 
tinuum shape. The model is successful in reproducing the 
soft X-ray excesses observed across a range of type-1 AGN 
(Middleton et al. 2007), together with their r.m.s. spectral 
variability (Gierlihski & Done 2006). The GD04 and SD06 
models are favoured by theoretical models of X-ray irradi- 
ated material in some sort of pressure balance, which show 
that the ionization instability limits the column density of 
material that can be present in the partially ionized state 
(Chevallier et al. 2006). This limited column is consistent 
with the optically thin absorption/emission models but is 
insufficient to dominate the entire photosphere of the ac- 
cretion disc as required for the reflection origin (Done & 
Nayakshin 2007). 

However, despite their success in modelling the data, 
the absorption/ emission models have some serious deficien- 
cies, particularly concerning the assumed velocity and den- 
sity structure of the wind, which must introduce consider- 
able scepticism in the validity of any conclusions based on 
these models (SD06). 

In this work we present a new approach to modelling 
the X-ray spectra of AGN observed through an accretion 
disk wind, that removes several of the limitations of previ- 
ous models. Specifically, we develop a code (hereafter XS- 
CORT; XSTAR Simulation Chain for Outflows with Radia- 
tive Transfer) that allows us to impose an external velocity 
field on a photoionization code in a self-consistent manner, 
and show that this results in large changes in the shape of the 
resulting X-ray spectrum. In Section 2 we present the details 
of the XSCORT code. In Section 3 we explore the behaviour 
of XSCORT spectra based on the average wind parameters 
from the sample of PG quasars presented in Middleton et 
al. 2007. We examine how the model spectra change as the 
outflow velocity field, total column density and initial ion- 
ization parameter of the wind are varied. In Section 4 we 
examine the impact that a complex ionization structure in 
the wind has on the XSCORT spectra and finally, in Sec- 
tion 5, we present our conclusions. 



2 XSCORT METHODOLOGY 

The GD04 and SD06 models of AGN X-ray spectra observed 
through an outflow rely on radiative transfer codes to calcu- 
late the absorption and emission signatures associated with 
partially ionized material. The results of these calculations 
are then convolved with an external velocity field to pro- 
duce a model X-ray spectrum that includes the impact of 
the accretion disk wind. The radiative transfer calculations 
are typically carried out under some simplifying assump- 
tions, the two most serious of which are that the gas is (a) a 
continuous, thin, slab of gas and (b) has a uniform density 
throughout the material. In addition to these, the externally 
applied velocity field necessarily represents a vast simplifica- 
tion of what is likely to be a considerably complex physical 
situation. In particular, the GD04 & SD06 models assume 
an external Gaussian velocity field, characterised by a very 
large dispersion (a ~ 0.2 — 0.3c). Such an external velocity 
field is only physically appropriate if the absorber consists 
of multiple, individual small clouds along the line-of-sight 
(l.o.s). In more continuous material, the velocity field sets 



the equivalent widths of the lines, so the photoionization 
cannot be calculated separately from the dynamics. 

In order to address these limitations we adopt a numeri- 
cal computation method that splits a single run of the radia- 
tive transfer code XSTAR into a series of multiple, linked, 
calculations that incorporates the effect of a global veloc- 
ity field on the calculations in a self-consistent way. XS- 
CORT also allows us to relax several of the other assumptions 
present in the previous models, such as constant density, giv- 
ing a considerably more realistic picture of the effects of an 
accretion disk wind on the X-ray spectra of AGN. 

2.1 Chaining XSTAR 

XSTAR calculates the absorption and emission spectrum 
that results from a stationary spherical shell of photoion- 
ized gas at a given radius from a given ionizing continuum, 
for a given a set of physical gas properties (density, abun- 
dance etc). XSTAR calculates the ionization balance and 
thermal equilibrium in a series of 'sub-shells' of gas, propa- 
gating the input incident spectrum through these sub-shells 
in turn, calculating the ionization balance, temperature and 
opacities in each. Absorption and emission lines are calcu- 
lated with an assumed Gaussian turbulent velocity broaden- 
ing added in quadrature to the self-consistently calculated 
thermal line widths. The turbulent velocity is normally the 
dominant contributor to the line widths, controlling the to- 
tal opacity in the absorption lines and hence the relation 
between the total column density of a given ion species and 
the observed line equivalent width (i.e. the curve of growth). 
It is important to note that the equivalent width of a line 
saturates as the line becomes optically thick, limiting its 
value, however a larger velocity width means that the in- 
trinsic line width is larger and the line is thus less likely to 
become optically thick at its core. As a result, a large turbu- 
lent velocity, for a given column, leads to larger equivalent 
width absorption line. 

XSTAR currently only implements this turbulent veloc- 
ity broadening on the absorption and emission lines, with the 
result that edges remain sharp despite the turbulent veloc- 
ity being characteristic of the material as a whole (SD06). In 
fact, even if this were calculated self-consistently, the Gaus- 
sian nature of the turbulent velocity field would still be a 
poor representation of the velocity field in a relativistic out- 
flow, where we would expect blue-shifted absorption along 
the l.o.s and both red and blue-shifted emission from all 
possible l.o.s. 

The computational method employed in XSCORT is 
similar to that used internally in XSTAR in that we divide 
the total column up into a large number of shells, each with 
smaller column densities, and treat each shell in a sepa- 
rate, albeit linked, manner. In addition to this, however, 
the effects of the outflow velocity field are incorporated self- 
consistently for each shell in the chain. Specifically, we step 
through the shells in order using the relativistically velocity 
shifted output from the XSTAR calculation in the (i-l) th - 
shell as the input for the XSTAR calculation in the i -shell, 
combined with a set of pre-defined gas parameters for the 
i -shell. Typically we specify five inputs for each shell in the 
XSCORT chain; column density (Nh,0i g as density (n;), ion- 
ization parameter (£»), radial velocity (v ra d,i) and the input 
ionizing spectrum. 
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XSTAR operates in the rest-frame of the absorbing ma- 
terial, so the incorporation of the outflow velocity requires 
that for each shell in the chain, the input luminosity, den- 
sity and radius all be transformed into the rest-frame of the 
shell, using the appropriate relativistic relations, in order for 
XSTAR to calculate the ion populations correctly (see Ap- 
pendix A2 & A3). In addition, there is a subtlety in the way 
XSTAR treats line and continuum escape probabilities which 
means that it is not appropriate to simply take the output 
luminosity from a single XSTAR run and use this as the input 
luminosity to the next calculation. Instead we use the same 
photon field XSTAR propagates internally (Appendix Al) so 
that XSCORT is consistent with a single, constant density, 
XSTAR run, in the absence of an external velocity field and 
constant gas density throughout the chain. 

The velocity shear within each shell in the chain is char- 
acterised using the turbulent velocity parameter in XSTAR, 
Vturb,i- We set the turbulent velocity in shell i, v tU rb,i to 
be the standard deviation for a uniform distribution, with 
an additional component, vt, representing the intrinsic real 
turbulent velocity of the material: 

Vturb,i = (Vrad,i — V ra d,i-l)/Vl2 + V t (1) 

For the spectra presented in Sections 3 & 4 the intrinsic 
turbulent velocity, vt, is set at 100 km s" 1 . The derivation of 
Vturt,i is given in detail in Appendix A4, but intuitively, it is 
clear that the Gaussian dispersion must be less than half the 
velocity shear because a Gaussian distribution is symmetric 
about the mean. 

2.2 Calculating the transmitted and emitted 
spectra 

The transmitted spectrum from shell i is calculated by trans- 
forming the transmitted spectrum from the i — 1 th shell, in 
the observers frame, into the rest-frame of shell i, absorbing 
the spectrum with the rest-frame opacities calculated by XS- 
TAR, and then transforming the newly absorbed spectrum 
back into the observers frame. Iterated over all the shells in 
the chain this then produces the total transmitted spectrum 
from the outflowing wind. 

Once in the rest-frame of shell i, the transmitted spec- 
trum is calculated in a similar fashion to the calculation XS- 
TAR performs internally (see Appendix A). At this point, 
however, we deviate from our requirement of consistency 
with a single XSTAR run by including a simple treatment 
of electron scattering (something XSTAR does not currently 
incorporate), scattering flux out of the transmitted (l.o.s) 
spectrum and into the emitted (global) spectrum. We follow 
the notation used in the XSTAR manual* and denote all spe- 
cific luminosities (units of erg/s/erg) as L t . The transmitted 
spectrum from shell i, denoted as L h e i in the XSTAR manual, 
is then given by: 

Lli = Ll i _ 1 e- (Tt - t - i+NH ^ T) (2) 

where L\ i _ 1 is the output transmitted spectrum from the 
previous shell, transformed into the rest-frame of shell i and 

* http://heasarc.gsfc.nasa.gov/docs/software/xstar/docs/ 
html / xstarmanual. html 



Ttot,i is the sum of the bound-bound, bound-free and free- 
free optical depths, tt = Nhjctt is the optical depth to elec- 
tron scattering, where Nh,% is the column density of the shell 
and or is the Thompson cross-section. We note that this is 
simplistic treatment of electron scattering, neglecting the 
effects of absorption, multiple scatterings, down-scattering 
and the Klein-Nishina cut-off. The latter two of these ef- 
fects are largely unimportant at energies <100 keV, while 
the contribution of multiple scatterings depends on Nhj but 
has less than a 10% effect on the spectral shape for column 
densities of ~10 23 cm" 2 . 

While the impact of the velocity on the transmitted 
spectrum is relatively simple to calculate, the impact on the 
emission spectrum is somewhat less straightforward due to 
the global nature of the wind material. We assume that the 
emission originates from a spherically symmetric wind, and 
thus the total emission spectrum will include a contribution 
from every possible l.o.s. We divide the spherical shell, i, into 
segments of solid angle dQ = sin 8 dOdcj>. The radial emission 
spectrum from each solid angle segment is then the emission 
spectrum output from XSTAR, re-normalised by dQ/4-K. Fi- 
nally, the total l.o.s emission spectrum from shell i is then 
the sum of the emission spectra from each solid angle seg- 
ment, shifted and beamed by the component of the outflow 
velocity of the segment along the l.o.s. The emission from 
each solid angle segment will be absorbed by wind material 
along its l.o.s to the observer (which will vary from almost 
nothing for material in the outer shells on the near side of the 
source, to much more than the continuum l.o.s absorption 
for material on the far side of the source) however, calculat- 
ing the absorption for the emission spectrum increases the 
computational complexity of the code enormously, and dra- 
matically increases the run-time of the program, so for the 
work presented here we restrict ourselves to calculating the 
unabsorbed total emitted spectrum and stress that this is 
an overestimate of the emission from the wind. 

The emission spectrum is also sensitive to both the lo- 
cal and global covering fraction of the wind material. Reduc- 
ing the local covering fraction used by XSTAR (i.e. whether 
the material in the immediate environment of each l.o.s is 
clumpy or uniformly smooth) can reduce the overall normal- 
isation of the emission spectrum by a factor ~2, but does 
not systematically change the shape of the emitted spectrum 
(we hardwire the local covering fraction C/j OC ai=0.001). Al- 
ternatively, imposing a large scale geometry on the wind ma- 
terial (evacuating a bi-conical structure along the spin axis 
of the central black hole, for example) can potentially have 
a considerable effect, but only on the emission spectrum. 
In the case of a bi-conical structure, a continuous region of 
solid angle is removed as a source of emission, resulting in a 
change in the detailed shape of individual emission features 
as well as a reduced overall normalisation. In the most ex- 
treme case, where the only remaining material is distributed 
in a thin equatorial disk, this results in the classic double- 
peaked emission lines familiar from accretion disk studies 
(e.g. Chen & Halpern 1989). The model spectra presented 
are calculated for a global covering fraction, Cf, g iobai=l, and 
thus represent the maximum possible contribution of the 
emission component to the total spectrum. 

The rest-frame emitted spectrum is calculated in a simi- 
lar fashion to the calculation XSTAR performs internally (see 
Appendix A), however we again deviate from our require- 
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ment of consistency with a single XSTAR run, by adding the 
electron scattered component to the genuine emission spec- 
trum. Again, using the same notation as the XSTAR manual, 
the total rest-frame emission spectrum from shell i, denoted 
Lj i, is thus given by: 

— L tt bb,i + L et bf,i + L t jf t i 

_i_ r, 5 . , P -{^H,i"T ) I , cr\ 

~ -t-'Ci—l I outward V / 

where L tl bb,i, L e ,bf,i & L e jf,i are specific luminosity spectra 
from bound-bound, bound-free and free-free atomic transi- 
tions. We note again that tt = NhjUt is the optical depth 
to electron scattering, where Nh,i is the column density of 
shell i, (tt is the Thompson cross-section and thus the final 
term in this equation represents the flux scattered into our 
l.o.s through electron scattering. 

The transmitted and emitted spectra are output on the 
same energy grid as that used by XSTAR, namely a loga- 
rithmic energy grid spanning E=0.1-10 6 . 

2.3 Ionization Structure 

The relativistic nature of the outflow has a strong impact on 
the ionization structure of the wind. For an outwardly ac- 
celerating wind, the input spectrum for shell i (L\ i _ 1 - see 
Appendix Al) will be both red-shifted and beamed (away) 
relative to shell i — 1. These velocity transformations (Ap- 
pendix A3) have a strong impact on the luminosity of the 
input ionizing spectrum (a factor of ~3 for terminal veloci- 
ties of ~0.3c), reducing the ionization state of the wind ma- 
terial at each step. The density of the wind material along 
the l.o.s is also affected by the wind outflow velocity. Specif- 
ically, the density denned in the observers frame is related 
to the density in the rest frame of the shell through a length 
contraction based on the l.o.s outflow velocity. The density 
is reduced by a factor of I/7, where 7 = (f — wf/c 2 ) 1 ^ 2 and 
Vi is the l.o.s component of v ra d,i, the radial velocity of shell 
i. The impact of the length contraction is relatively small 
(~5% for v out — 0.3c) in contrast to the impact that the 
transformed luminosity has on the ion populations. 

Additionally, imposing an external density structure on 
the wind will impact on the ionization structure of the wind. 
In order to understand the properties of the XSCORT mod- 
els, in Section 3 we impose a density structure on the wind 
material such that the changes in density counteract the 
red-shift and beaming effects that affect the input ionizing 
spectrum. To first order this results in a constant ionization 
state throughout the material; a situation that is broadly 
familiar from static models of ionized material. The remain- 
ing change in the ionization parameter (typically <3%) as 
we move through the wind material reflects only the de- 
crease in luminosity due to absorption by the wind material 
and the increase in distance of each shell from the ionizing 
source, due to the geometric width of the wind. 

The ionization instability provides a natural physical 
mechanism for producing a wind with a constant ionization 
state structure. Moving through a slab of highly ionized gas, 
the transition from highly ionized material to low ionization 
material occurs over a small range of radii. As the ionization 
state of the material drops, however, the material will cool 
and collapses, increasing its density which, in turn, drives 
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Figure 1. Top Panel: The UV/X-ray spectra (E Le) that re- 
sults from an XSCORT run with a linearly accelerating veloc- 
ity field. Bottom Panel: The X-ray spectra that results from an 
XSCORT run without an external velocity field. Upper: The to- 
tal spectrum (dark green) and the input power-law continuum 
(black). Lower: The transmitted spectrum (blue), the emitted 
spectrum (red) and the input power-law continuum (black). 

the ionization state down even further. This runaway pro- 
cess results in almost neutral, dense, clumps of material, 
however their volume filling fraction is small and as a result 
these clumps are not expected to have a large impact on 
the X-ray spectrum. The dominant contribution to the X- 
ray spectra shape is, then, from the relatively thin layer of 
partially ionized material with (almost) constant ionization 
state (Chevallier et al. 2006). 

This scenario does not, however, account for the im- 
pact of the outflow velocity. The clumps will collapse on 
timescales of approximately sound crossing time of the 
clump. If the wind outflow velocity is considerably greater 
than the sound speed in the gas before it begins to collapse, 
then these clouds may not be able to collapse effectively. 
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Figure 2. Velocity as a function of step through the XS- 
CORT chain. The red line corresponds to linear acceleration 
(p=l). The dark and light green lines have p=0.2 & 0.5, respec- 
tively and represent the situation where the majority of the accel- 
eration of the wind happens near the launch radius. Conversely, 
dark and light blue lines have p=2 & 5 and represent the situa- 
tion where the majority of the acceleration of the wind happens 
far from the launch radius, in each case, the launch velocity is 
and the outer-edge velocity is 0.3c. 

In Section 4 we relax the previous constraint on the den- 
sity structure of the wind and explore winds with a density 
structure that thus have wide range of ionization parame- 
ters through the wind material. We take the simplest possi- 
ble approach to the density structure of the wind and treat 
only structure in the radial direction. Because the transmit- 
ted spectrum only samples a pencil-beam l.o.s, it is only 
sensitive to the radial density structure of the wind and is 
insensitive to the covering fraction or the filling fraction of 
the gas. In contrast, the emitted spectrum will be strongly 
affected by the 3D density structure of the flow, however the 
emission is not the dominant component of the total spec- 
trum except in calculations with particularly high column 
densities, or low ionization states. Calculating the impact of 
a detailed 3D structure on the emitted spectrum is beyond 
the scope of the work presented here and, for the majority 
of the spectra presented here, is not expected to strongly 
impact the shape of the resulting spectra. We stress again 
that the spectra presented here have the maximum possible 
contribution of the emission component to the total spec- 
trum (see Section 2.2) and as such should be treated as an 
upper limit in this regard. 



3 BEHAVIOUR OF THE OUTPUT SPECTRA 
FROM XSCORT 

In this section we examine the behaviour of the XS- 
CORT models as key parameters are varied, to highlight the 
impact that specific parameters have on the overall shape of 
the X-ray spectrum seen through an accelerating wind. We 



use a power law continuum as the initial ionizing spectrum 
and we define a standard set of parameters based, in part, 
on the set of 'average' AGN wind properties from Middleton 
et al. 2007. 

Middleton et al. (2007) analyse X-ray observations of 
a moderate sized sample of nearby, bright, PC quasars 
and NLSls, and show that the X-ray spectra of this sam- 
ple can be equally well fit by partially ionized, optically 
thin, smeared, material seen in absorption, or by optically 
thick, partially ionized, material seen in reflection (Crummy 
et al. 2006). The average absorption in the Middleton et 
al. sample has an ionization of logio(£)=3, a column den- 
sity of N_sr=3xl0 23 cm -2 and an outflow velocity of >0.3c, 
although there is a considerable range across the sample in 
each of these parameters (factors of >10 in the first two 
cases and a factor of ~ three in the latter cases) . The stan- 
dard set of parameters is then: a 100 step chain, Logio(£)=3, 
N^ iiot =3x 10 23 cm -2 , an initial gas density of 10 12 cm" 3 
and linear acceleration between and 0.3c. In each of the 
following subsections we then vary one of these parame- 
ters. With these parameters, the radial depth of the winds 
simulated with these parameters is small compared to the 
launch radius. However the physics that governs the shape 
of the output X-ray spectra depends only the distribution 
of column density and atomic species as a function of out- 
flow velocity. Under the assumption that the gas is smooth, 
reducing the gas density in order to simulate a geometri- 
cally thicker wind, would thus force the launch radius of 
the wind out to large radii in order to have the given ion- 
ization parameter but the resulting spectrum would stay 
approximately the same. These parameters are also similar 
to those resulting from hydrodynamic simulations of accre- 
tion disc winds in AGN, though these also have a complex, 
filamentary structure which is beyond the scope of this pa- 
per. Nonetheless, these filaments have gas densities ~10 9-13 
cm -3 and maximum velocities of ~ 0.2c, with the major- 
ity of the acceleration occurring close to the base of these 
structures (Proga, Stone & Kallman 2000; Proga & Kallman 
2004). 

3.1 The outflow velocity field 

We examine the effect of a basic accelerating velocity law 
on the resulting X-ray spectrum, adopting a simple p-type 
form; i.e. the radial velocity of shell i is given by: 

Vrad,i = Vl + (V - Vl) ( -r^— ) (4) 

where vi and v are the launch and terminal velocities, 
respectively, and i m ax is the number of steps in the XS- 
CORT chain. We note that for p=l, equation 4 produces a 
linear acceleration through the wind material. 

Figure 1 shows an example of the X-ray spectrum that 
results from an XSCORT run simulating an accelerating 
wind, over the 7 eV to 13 keV range, and a spectrum that 
does not include an external velocity field. The contrast be- 
tween the two spectra highlights the massive impact that 
the velocity field has on both the transmitted and emitted 
spectrum. The XSCORT spectrum is clearly dominated by 
two prominent absorption regions around ~0.06 and ~1.5 
keV. The first of these absorption conglomerates is the re- 
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Figure 3. The 0.3-13 keV X-ray spectra (E Le) that result from XSCORT runs with different p velocity laws: Panels 1-3 (Top, 
left to right): p = 0.2, 0.5 & 1.0 (linear acceleration). Panels 4-6 (Bottom, left to right): p = 2.0, 5.0 &; the Swindabs+Swindemm 
model from Schurch &; Done (2006). All other parameters are fixed at the standard values. Upper: The total spectrum (dark 
green) and the input power law continuum (black). Lower: The transmitted spectrum (blue), the emitted spectrum (red) and the 
input power-law continuum (black). We note that Middleton et al. (2007) fit only smeared absorption to their sample, ignoring 
the associated emission component. Thus, the blue line in the bottom right panel represents the smeared absorption Middleton 
et al. fit. We stress that absorption of the emission spectrum in subsequent shells is neglected throughout this work and that all 
the spectra are calculated for a global covering fraction of Cf g i i, a i=l. Thus, the emission spectra presented here represent the 
maximum possible contribution from emission to the total spectrum. 



suit of low ionization species (CIII, CIV, NIII, FeV etc) and 
falls well outside the band-pass of current X-ray observato- 
ries. Interestingly, the UV absorption is similar to the broad 
absorption lines that are characteristic of BAL QSO's, a 
possibility that we will examine in future work by combin- 
ing XSCORT with hydrodynamic simulations of accretion 
disk winds. The second of the absorption troughs is the re- 
sult of higher ionization species (NVI, OVII, OVIII, FeXVII, 
etc) and occurs in a region that is well covered with the cur- 
rent generation of X-ray observatories. While the stationary 
model has the continuum absorption dominated by bound- 
free edges, with individual absorption lines superimposed on 
this, the increased velocity shear in the wind model results 
in the resonance lines completely dominating the opacity. 
This produces a dramatic increase in the apparent absorp- 
tion of the continuum spectrum when a velocity field is in- 
cluded in the calculations. For the remainder of this work 
we concentrate of the 0.3-13 keV band-pass observable with 
XMM-Newton & Chandra, although we note that the ~0.02- 



0.09 keV absorption feature may well be detectable in cur- 
rent, or future, far-UV spectra of AGN. 

Figure 2 shows the wind velocity as a function of frac- 
tion through the total processing chain, for several values of 
p. There is a considerable range in the mean velocity of the 
simulated winds in Figure 2 (< v >=ii o /(p+l)=0.05 & 0.25c 
for p=5 & 0.2) however the velocity dispersion is similar in 
each case (<r=< v > p/(2p + l) 1/2 =0.07c & 0.04c for both 
p—5 & 0.2). Figure 3 (panels 1-5) shows the model X-ray 
spectra from XSCORT with these velocity laws, highlighting 
the range of spectral shapes that result. To zero-th order 
the transmitted component of the X-ray spectrum remains 
largely insensitive to the details of the velocity field because 
this component is dominated by the effects of velocity dis- 
persion. In contrast, the emission spectrum, which includes 
contributions from red-shifted material from the far side of 
the wind as well as blue-shifted material from the near side 
of the wind, is sensitive to a combination of both the mean 
velocity and velocity dispersion. 
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Figure 4. Structure in the iron K« region as a function of velocity 
law. The dark green, green, red, dark blue & light blue lines corre- 
spond to the velocity laws shown in Figure 2 (p=0.2, 0.5, 1.0, 2.0 
& 5.0, respectively). The black line is the input power-law contin- 
uum. Upper Panel: Emission spectra. Lower Panel: Transmission 
spectra. High p- values result in narrow iron line and iron edge fea- 
tures at almost rest-frame energies. Low p- values, however, result 
in considerable smearing, particularly of the iron line, producing 
features with maximum emission and absorption at energies that 
correspond to the maximum velocity of the outflowing material. 

Detailed differences in the shape of isolated emission, 
and absorption, features such as the iron Ka line and edge 
(Figure 4) explicitly reveal the distribution of the column 
density of a specific ion with velocity. At high p-values (e.g. 
p=5) a considerable amount of the wind material is present 
in regions where the velocities and, crucially, the velocity 
gradient are both small. For absorption features, sensitive 
only to the velocity dispersion, most of the opacity at high 
p-values occurs at near-rest-frame energies with an tail to 
higher energies. The emission features from a wind with a 
high p-value are not highly smeared and are imprinted on 
the spectrum as narrow features with maximum absorption 
and emission at almost rest-frame energies. Conversely, at 
low p-values (e.g. p=0.2) a considerable amount of the wind 
material is present in regions where the velocities and, cru- 
cially, the velocity gradient are both large. For absorption 
features, this results in the majority of the opacity occurring 
at energies corresponding to the mean blue-shift of the ma- 
terial, with a tail to lower energies. For the emission features 
this results in strongly broadened features with maximum 
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Figure 5. The ratio of the standard linearly accelerating XS- 
CORT model, convolved with the XMM-Newton EPIC PN in- 
strument response, to a simple power-law, also convolved with 
the instrument response (green). Also shown is the ratio of the 
smeared ionized absorption used in Middleton et al. (2007) con- 
volved with the same instrument response, to a simple power-law 
(blue). 

emission at energies that, again, correspond to the mean 
blue-shift of the material (Figure 4, upper panel). 

This highlights an issue concerning the strong iron Ka 
features observed in some nearby Quasars. These model 
demonstrate that absorption in an accretion disk wind can 
make the sharp drop seen around 7 keV in some AGN, as 
long as the majority of H-like iron column is at low velocity 
(Done et al. 2006). However, the emission from this material 
is then far too narrow to create a smooth P-Cygni profile 
such as we might expect from a wind. The impact of the 
emission can be reduced by reducing the covering factor of 
the wind, in which case the profile is then dominated by the 
absorption, but we caution that our emission model does 
not include resonant line scattering which will enhance the 
line emission in the emitted spectrum. A full treatment of 
scattering, absorption and emission of resonance lines is ex- 
tremely complex and beyond the scope of this work (Matt 
et al. 1996). 

More importantly, these models show a key difference 
from both the observed shape of the soft X-ray excess and 
the simplistic, Gaussian-smeared wind models. Firstly, for 
the same column density of material, all the accelerating 
wind spectra show considerably more absorption than the 
spectra from a slab ionization model convolved with a Gaus- 
sian velocity field. The proper treatment of velocity in these 
models means that the absorption lines stay on the opti- 
cally thin part of the curve of growth and their equivalent 
width, for the same column, is thus considerably greater. 
We note that it is not quite equivalent to simply use a lower 
column in the XSCORT models, because the impact of the 
wind velocity enhances the line absorption over the bound- 
free edges. 

Even more importantly, the accelerating wind spec- 
tra have sharp features in the total spectrum that are not 
present in the GD04 & SD06 models. The presence of narrow 
features is a product of the reasonably small velocity dis- 
persion in the velocity laws shown, compared to the ~0.3c 
velocity dispersion required by the GD04 and SD06 mod- 
els. The sharp features predicted by the XSCORT models 
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Figure 6. As for Figure 3 but for XSCORT runs with different total column densities: Panels 1-3 (Top, left to right): Njj 



9xl0 21 , 3xl0 22 & 6xl0 22 cm- 
are fixed at the standard values. 



. Panels 4-6 (Bottom, left to right): N ff = 9xl0 22 , 3xl0 23 , 6xl0 23 cm -2 . All other parameters 



are observable even with current high resolution X-ray CCD 
detectors. Figure 5 shows the total spectrum from the stan- 
dard linearly accelerating XSCORT model convolved with 
the XMM-Newton EPIC PN instrument response divided 
by a simple power-law model, also convolved with the in- 
strument response. Also shown is the ratio for the smeared 
ionized absorption used in Middleton et al. (2007). The over- 
all shape of the smeared ionized absorber model and the XS- 
CORT wind models are broadly similar, albeit that the XS- 
CORT wind model shows greater absorption, for the same 
column density. However, the sharp features at 0.6 & 1.4 
keV are both clearly visible in the XSCORT model, as is the 
strong drop from H-like iron Ka absorption lines around 7 
keV. The strong, line-like, feature at ~1.4 keV is, in fact, 
not a line (although there is a considerable contribution to 
the flux at these energies from line emission). It is primarily 
the result of considerable broadened absorption at energies 
directly above and below 1.4 keV, associated with ionized 
oxygen and iron respectively. While it is currently not pos- 
sible to fit the XSCORT models to observed X-ray spectra, 
the presence of this dramatic feature suggests that the XS- 
CORT models presented here will be unable to reproduce 
the observed smooth featureless soft X-ray excess. Again, 
we stress that the effect of absorption on the emission spec- 
trum has been neglected in the calculations presented here, 
and that this will imprint additional features on the emission 



spectrum and thus have a considerable impact on the total 
spectrum. However we note that these features will be simi- 
lar in character to the features already evident in the trans- 
mitted spectrum and, as such, will not remove the sharp 
features that are present in the current total spectrum. In 
fact, in regions of the spectrum where the (comparatively) 
smooth emitted spectrum dominates the total spectrum, in- 
corporating the absorption on the emission spectrum will 
introduce additional sharp features into the resulting total 
spectrum! Clearly this will not result the XSCORT models 
being more compatible with the simplistic smeared ionized 
absorption models or the observed X-ray spectra of AGN. 

The sharp features, particularly the feature at ~1.4 
keV, can be removed, however, by increasing the terminal 
velocity of the wind, increasing the energy range over which 
the prominent absorption troughs have a strong impact. The 
mean velocity dispersion required, ~0.3c, corresponds to a 
terminal velocity of ~c, implying an association with ma- 
terial in the jet rather than material in an accretion disk 
wind. 

3.2 Total column density 

Figure 6 shows the X-ray spectra from winds with a range 
of total column densities. 

At low column densities (<6x 10 22 cm -2 ) the total spec- 
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trum is dominated by the transmitted spectral component. 
Above these columns the absorption is significant enough 
that the emission spectrum begins to dominate regions of the 
spectrum. At very high column densities (>fewx 10 23 cm -2 ) 
the absorption is so significant that the <10 keV spectrum 
becomes completely dominated by the emission component. 
We stress, again, that the emission spectrum calculated by 
XSCORT is, currently, not self-consistently absorbed by the 
wind material, resulting in an emission spectrum that is 
somewhat brighter than we would otherwise expect, partic- 
ularly at high column densities (see Section 2.2). For most of 
the spectra shown here, the ionization parameter does not 
change appreciably across the thin wind currently simulated 
by the XSCORT chain. As such, the ionic species present in 
the wind do not change considerably as a function of ra- 
dius and, thus, the absorption and emission features do not 
change character significantly, despite increasing consider- 
ably in strength. For thicker winds (or columns >fewxl0 23 
cm -2 ) the ionization parameter will change considerably 
across the width of the wind and the resulting changes in 
the ionic species with radius results in considerable changes 
in the shape of both the emitted and transmitted spectrum 
(see Section 4). 



3.3 Ionization parameter 

Figure 7 shows examples of the X-ray spectra from winds 
with a range of initial ionization states. 

At low ionization states (Logm(£)<2.5) the X-ray spec- 
trum is dominated at high energies by the transmitted 
spectrum and at low energies by the emission spectrum 
(again, somewhat overestimated by lack of self absorption). 
Above (Logio(£)>2.5) the total spectrum becomes increas- 
ingly dominated by the transmitted component as there is 
less absorption and hence less emission, and the features im- 
printed on the X-ray spectrum become dominated by high 
ionization species. In particular, the iron Ka absorption line 
becomes stronger as more of the column is ionized to He- 
and H-like iron. 



4 A COMPOSITE COLUMN 

In the winds we have simulated so far, the ionization parame- 
ter does not change significantly through the wind material 
unless the column densities are extremely large. However 
the gas properties (most notably the density structure - see 
Section 3.2) in these calculations may not be an accurate 
description of the real physical conditions in an AGN wind. 
In the more case of a geometrically thick wind, where ei- 
ther the turbulent velocity, or the outflow velocity shear, is 
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Figure 8. Density as a function of step through the XS- 
CORT chain. The black line corresponds to constant density; the 
red line corresponds to linearly decreasing density. The green and 
blue lines are q-law profiles, similar in form to those used for the 
velocities in Figure 2. The green lines correspond to g=0.2 & 0.5 
(dark and light, respectively) and the blue lines correspond to q=2 
& 5 (again, dark and light, respectively). In each case, the initial 
density is 10 12 cm -3 and, in the non-constant density cases the 
outer-edge density is 10 10 cm -3 . 



greater than the sound speed in the gas, a wide range of 
ionization states can be present throughout the material. 

The variation in the ionization parameter through such 
a composite column depends on the velocity and density 
structure of the wind (Section 2.3) and on the increasing 
distance of the gas from the source of the ionizing radia- 
tion. We investigate the spectra from a range of composite 
columns by imposing a range of density profiles on the wind 
material. Along with a simple constant density profile, we 
examine density profiles similar in form to the velocity pro- 
files presented earlier. With the exception of the constant 
density case, the density of shell i in the rest frame of the 
observer is given by: 



Figure 9, panel 1, with the X-ray spectrum from the stan- 
dard set of wind parameters (Figure 3, panel 3), clearly high- 
lights the impact that composite columns have on the X-ray 
spectrum. In the constant density case, the ionization pa- 
rameter drops by a factor of ~4 across the wind due to 
the reduction in luminosity of the source due to relativistic 
beaming which is now not matched by a change in density to 
keep the ionization parameter constant. Thus the constant 
density wind includes material at lower ionization state, giv- 
ing considerably more absorption on the transmitted spec- 
trum than in the standard case. The change in the ion pop- 
ulations is particularly evident in the shape of the iron K 
absorption line, where the absence of high ionization species 
at high velocity results in a more 'edge-like' shape compared 
to the rectangular notch in the standard spectrum. 

Panels 2-6 show the changes in spectral shape that re- 
sult from strongly decreasing density profiles with different 
q- values. In these cases the ionization structure is more com- 
plex than in either the standard spectrum case (where the 
ionization state remains roughly constant throughout the 
wind), or the constant density case. In the two high-g cases, 
the density remains high for a considerably amount of the 
total column. While the density remains high the ioniza- 
tion parameter drops (in a similar fashion to the constant 
density case), however once the density drops considerably 
(O~20%) the ionization state of the material becomes dom- 
inated by the density of the material and the ionization pa- 
rameter begins to rise. The final ionization state at the outer 
edge of the cloud is higher than the initial ionization state. 
In the most extreme high-g case shown, the ionization state 
rises by a factor of ~4.5. In the low-g cases, the density drops 
rapidly and ionization state of the material is dominated by 
the dropping density throughout the wind, rising rapidly. In 
the most extreme \ow-q case shown the ionization state at 
the far edge of the wind is ~24 times higher than the initial 
ionization state. 

The spectra from composite columns are, in fact, very 
similar to the spectra calculated with (approximately) uni- 
form ionization state (compare Figure 9 with Figure 7), 
demonstrating that complex ionization structure columns 
can be modelled with more simple models that have uni- 
form ionization state (at least to first order). The best-fit 
ionization state from these fits represents an average for the 
ionized material. 



Hi = ni + (n - n ; ) I J (5) 

where n; and n are the launch and terminal densities, 
respectively, and imax is the number of steps in the XS- 
CORT chain. We note that, in a similar vein to the velocity 
profiles presented earlier, there is no specific physical moti- 
vation for any of the density profiles used here. Rather, the 
resulting spectra highlight the general impact that a com- 
posite column has on the absorbed X-ray spectrum of AGN. 

Along with the constant density case, we calculate 
model spectra for five g-law density profiles (q = 5.0, 2.0, 
1.0, 0.5, 0.2), shown in Figure 8. Each profile, with the ex- 
ception of the constant density case, has a launch density of 
10 12 cm -3 and a terminal density of 10 10 cm -3 . 

Figure 9 (panels 1-6) shows the model X-ray spectra 
from XSCORT with these density profiles. A comparison of 



5 CONCLUSIONS 

We adopt a numerical computational method that links a se- 
ries of radiative transfer calculations, including the effect of 
a global velocity field in a self-consistent manner (XSCORT), 
in order to produce a considerably more realistic picture of 
the X-ray spectrum of an AGN seen through an accretion 
disk wind. 

We present a series of example spectra from the XS- 
CORT code that allow us to examine the shape of AGN 
X-ray spectra seen through a wind for a range of veloc- 
ity and density distributions, total column densities and 
initial ionization parameters. The detailed spectral models 
presented here clearly show considerable complexity and de- 
tailed structure that is strongly affected by all these factors. 
Nowhere is the impact of these key parameters more clearly 
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Figure 9. As for Figure 3 but for XSCORT runs with different density profiles: Panels 1-3 (Top, left to right): constant density, q 
= 5.0 &; 2.0. Panels 4-6 (Bottom, left to right): q = 1.0, 0.5 & 0.2. All other parameters are fixed at the standard values. 



highlighted than in the region around iron Ka. In partic- 
ular, we use this region to highlight the changes in shape 
of the absorption and emission features as a function of the 
assumed velocity law (see Figure 4). 

The complexity and structure of the XSCORT models 
is not present in the ad-hoc spectral models, presented in 
GD04 & SD06. The presence of sharp features in the XS- 
CORT spectra contrasts strongly with the previous models 
and with the smooth nature of the observed X-ray spectra of 
AGN with soft X-ray excesses. We note, however, that the 
models presented here have an energy resolution consider- 
ably in excess of that of typical CCD detectors available on 
current X-ray satellites and that the sharp features in these 
models will be strongly smeared when convolved with an ap- 
propriate instrument response (see Figure 5). Despite this, 
the XSCORT models clearly show considerable differences 
in shape in comparison with the previous smeared ionized 
absorption models. In particular, the XSCORT spectra have 
a strong, line-like, feature at ~1.4 keV. This feature is, in 
fact, not a line (although there is a considerable contribu- 
tion to the flux at these energies from line emission), it is 
primarily the result of considerable absorption at energies 
directly above and below 1.4 keV, associated with ionized 
oxygen and iron respectively. 

Although the current version of the XSCORT code has 
some unaddressed issues, it is extremely unlikely that the 
presence of these strong features will be dramatically af- 



fected by addressing these points. The most significant pa- 
rameter affecting the presence of these features is the termi- 
nal velocity that the wind reaches. Increasing the terminal 
velocity, and hence increasing the velocity dispersion across 
the wind, will extend the region of the spectrum that is 
directly impacted upon by the ionized oxygen and iron ab- 
sorption. Thus, a terminal velocity of ~c could, potentially, 
result in a spectrum somewhat closer to the ad-hoc smeared, 
ionized absorption models, presented in GD04 & SD06, that 
are known to fit the observed data well. Such fast moving 
material cannot be associated with a radiatively driven ac- 
cretion disk wind, however such high velocity material may 
exist in a radio jet. 
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APPENDIX A: XSTAR DETAILS 

The inner workings of XSTAR have been presented else- 
where, notably in the XSTAR manual*, however here we 
reiterate several of the more subtle and important details 
that effect the XSCORT code. 



The most significant issue concerning XSTARs internal pro- 
cesses is the spectrum XSTAR propagates internally. The 
output transmitted and emitted, spectra from a single XS- 
TAR run (denoted , & L^ 7 ' , respectively, in the XS- 
TAR manual [chapter 8, section E.4]) are not the same quan- 
tity that XSTAR propagates between sub-shells in order to 
calculate the ionization balance and thermal equilibrium. 
Using the same notation as the XSTAR manual, the local ra- 
diation field propagated from sub-shell i into sub-shell i + 1 
is defined as: 



7" 1 T 1 



-(T*/,i+T-//, 4 ) 



(Al) 



where L 1 j is the specific luminosity output from the sub- 
shell, i (and hence the input to shell i + 1) h\ is the 
specific luminosity output from the previous, i — 1 sub-shell 
(and thus is the input to shell i - in the case of the initial 
shell, i=l, this is equal to the input incident spectrum, £ e ,o), 
Tbf.i & T ff,i are the optical depths of bound-free and free- 
free interactions in the sub-shell i, respectively, and £ e ,6/,i & 
i e ,//,i are, correspondingly, the specific luminosities emitted 
from sub-shell i due to bound-free and free-free interactions. 
By contrast, the final output spectra from a completed XS- 
TAR run are given by: 



,o e 



-(E 4 (Ti, / , 4 +T // , i +T I>M )) 



= £(*« 



+ u 



bf,i 



+ L 



(A2) 
(A3) 



L[ = £ {L €: bb,i + L e ^ b f,i + L e Jf ti ) j ou t W ard (A4) 
i 

where r bbti if the optical depth of bound-bound transitions 
in sub-shell i. Since we require consistency with a single 
XSTAR run, we must propagate the same thing that XS- 
TAR propagates internally, between our linked XSTAR runs. 
Unfortunately, it is not simple to reconstruct the locally 
propagated radiation field from the final sub-shell, i max , of 
an XSTAR run (L* ,; TOax ) using the output products of that 
XSTAR run (equations 2, 3 & 4). Instead, we modify the 
XSTAR source code directly to output i?iiox for a given 
XSTAR run, allowing us to use this quantity as the input 
spectrum for the next XSTAR run in the chain. 



A2 Ionization parameter 

Internally, XSTAR calculates the ionization state for each 
sub-shell based on the propagated specific luminosity from 
the previous sub-shell, the density of the sub-shell and the 
thickness of the sub-shell. Although possible in principle, we 
do not specify the radial profile of the ionization parameter 
for the XSCORT spectra presented here. Instead, we follow a 
similar process to the internal XSTAR calculation. The ini- 
tial ionization parameter is specified at the innermost edge 
of the wind material and then the ionization parameter in- 
put to each link in the chain is calculated as a function of 
distance and absorption along the l.o.s for each shell, i: 



* http://hcasarc.gsfc.nasa.gov/docs/software/xstar/docs/html/ 
xstarmanual.html 



(A5) 
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where Lj is the luminosity propagated into shell i (i.e. the 
luminosity of the specific luminosity spectrum propagated 
into shell i, L* calculated over the f-fOOO Rydberg range 
- see Equation Af), m is the rest frame gas density in shell i 
and Ri is the radius of shell i from the central X-ray source. 
Ri is given by: 

Ri = *+g(^) (A6) 

where Ri is the innermost launch radius of the wind, Nhj 
is the column density of shell j and nj is the density of the 
wind material in shell j. 

The radial dependence of the ionization parameter is 
largely unimportant here since the thickness of each shell 
in the winds simulated is much less than the launch radius 
of the wind, as is the thickness of the wind as a whole. The 
radial dependence of the ionization parameter will, however, 
be important for physical scenarios where the wind may have 
considerable radial thickness. 

A3 Velocity Transforms 

Velocity transforms are used calculate both the spectrum 
propagated between each shell in the chain, and the trans- 
mitted and emitted spectrum in the observers frame. In the 
case of a redshift the transformed specific luminosity spec- 
trum is given by: 

L' e = L^ 3 (l + v/cf (A7) 

where L e is the original specific luminosity spectrum, 7 = 
(1 — v 2 /c 2 ) 1 / 2 and v is the relative radial velocity of the two 
reference frames. 

In the case of a blue-shift, the transformed specific lu- 
minosity spectrum is given by: 

« - l -(w^m)' ,A8) 

A4 Velocity shear characterization 

We assume that the velocity shear across each individual 
shell is linear and we use the standard deviation of this dis- 
tribution as a representation of the actual velocty shear. The 
mean velocity in shell i is given by: 

< V ra d,i > = ^(V ra d,i-1 + Vrad,i) (A9) 

The mean squared velocity is given by integrating the 
velocity distribution as a function of column density through 
the shell (we note that for the following equations we have 
removed the subscript 'rod' on the velocity terms for clarity 
purposes): 



where Nn,i is the total column density of the shell i. The 
standard deviation is given by the standard equation of =< 
v rad,i > — < v ra d,i > 2 and substituting Equations A9 & 
All into this expression thus gives: 

<Ti = (V ra d,i - V ra d,i-l) / (A12) 

The turbulent velocity in the XSTAR run for shell i, 
Vturb,i IS then set to be equal to this standard deviation, but 
with an additional component, v t , representing the intrinsic, 
real, turbulent velocity of the material (see Equation 1). 



APPENDIX B: NUMERICAL LIMITATIONS 

Here we note some of the technical limitations of the XS- 
CORT simulation process. The number of steps in the chain 
and interpolation routines used for each of the velocity 
shifts, along with, the number of area segments in the emis- 
sion spectrum calculation, all present numerical limitations 
in the accuracy of the output spectrum. 

The most significant of these concern the number of 
steps in the chain and the interpolation errors associated 
with the velocity shifts. When the number of steps in the 
chain is small (<40 runs), the velocity grid is coarse and 
sinusoidal wiggles are introduced to the spectrum, particu- 
larly around discrete features. Conversely, when the number 
of steps in the chain is high (>150 runs), the cumulative 
errors from the large number of interpolations run over the 
course of the chain (three for each step) introduce significant 
changes in spectral shape of the output spectrum. In par- 
ticular, when the width of narrow features in the spectrum 
is comparable to the size of the energy grid, a large number 
of interpolations fills in the absorption features, and reduces 
the strength of emission features. In extreme cases, these 
errors can be very significant (factors of ~3-4) but they can 
be minimised either by ensuring that the turbulent veloci- 
ties used throughout the chain are always large enough that 
none of the spectral features are narrow in comparison with 
the energy grid, or by increasing the density of the energy 
grid prior to interpolation. The latter solution will be im- 
plemented in future versions of the code. 

By contrast, the number of area segments used in the 
emission spectrum calculation introduces relatively minor 
errors to the calculated spectrum. Using a small number of 
area segments (<200) results in discrete wiggles in the emis- 
sion spectrum, particularly on the blue-ward side of strong 
emission features. Using a large number of area segments 
(>5000) rapidly becomes prohibitively computationally in- 
tensive. 



<v~ > = 
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